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(N 

^ A two phase elastic composite with weakly compressible elastic inclusions is considered. The 

^T^ homogenised two-scale limit problem is found, via a version of the method of two-scale con- 

■^T vcrgence, and analysed. The microscopic part of the two-scale limit is found to solve a Stokes 

type problem and shown to have no microscopic oscillations when the composite is subjected 

>L/ to body forces that are microscopically irrotational. The composites spectrum is analysed and 

Cm shown to converge, in an appropriate sense, to the spectrum of the two-scale limit problem. 

A characterisation of the two-scale limit spectrum is given in terms of the limit macroscopic 

1 I and microscopic behaviours. 

T^ 1. Introduction 

H It is well known that composite materials often display physical properties that 

'— ' are not observed by their individual constitutive parts. This leads to the question 

,__( whether one can produce composite materials to exhibit a desired property. Stated 

^ differently, can one determine the effective properties of a composite material with 

^^ prescribed microscopic data. Mathematically, one can approach this question via 

the application of homogenisation theory to determine the 'homogenised' limit 



o 

f — . to the equations modelling the composite' s behaviour. Such limit equations can 



^ 



be then considered to contain the effective physical properties of the composite 
^^ with the limit solutions describing the effective behaviour. Historically, "Classi- 

cs cal" homogenisation has been applied to composites with moderately contrasting 

. . heterogeneity. Here the heterogeneity is replaced by an equivalent homogeneous 

K^ medium with uniform physical properties, and is therefore incapable of describing 

K^ a range of interesting and unusual effects. Subsequently, homogenisation was used 

?-H to study composite materials with highly contrasting coefficients. The so-called 

high contrast homogenisation theory has been used to describe many non-trivial 
and interesting behaviours; examples include memory effects (e.g. [1-3]) and other 
non-local effects (e.g. [4-7]). 

A useful analytical tool in the homogenisation theory is the method of two- 
scale convergence first introduced by Nguesteng [8] and substantially developed 
by Allaire [9] particularly in the context of high contrast periodic problems. Two- 
scale convergence was further developed by Zhikov for the study of high-contrast 
spectral problems in bounded [10] and unbounded [11] domains. Therein, Zhikov 
described a "two-scale limit operator" and the Hausdorff convergence of spectra 
in terms of strong two-scale resolvent convergence and the two-scale compactness 
of eigenfunctions. Zhikov also explicitly described the "limiting" spectrum by a 
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coupled system of limit equations in terms of the macroscopic and microscopic 
variables. Furthermore, Zhikov showed upon decoupling, the effective macroscopic 
properties to depend non- linearly on the spectral parameter, essentially giving rise 
to a description of a "microresonace" effect: a distinct change in physical properties 
when the applied 'macroscopic' frequency is close to the eigenfrequencies of the 
microscopic inclusions. In the context of electromagnetism these results first due 
to Zhikov, [10], where interpreted as the appearance of effective negative magnetism 
for appropriately polarised waves of certain frequencies in a non-magnetic material 
with high contrasting electric permeability, see [12]. 

Recently, the interest in high-contrast homogenisation has been boosted by ap- 
plications to so-called metamaterials, c.f. [13]. This class of composite materials ex- 
hibit macroscopic or 'effective' physical properties not commonly found in nature, 
such as electromagnetic meta-materials with negative refractive index, cloaking de- 
vices and super lens, to name but a few. More recently, Smyshlyaev following [4], 
showed in [14], using multi-scale asymptotic expansions, that composite materials 
with a "partially high contrast" between constitutive parts are capable of account- 
ing for additional phenomena such as directional propagation. Thus one may find 
in the study of this broader class of composite materials the physical effects one 
seeks in meta-materials. In [15] the two-scale homogenisation analytic tools are 
developed for a broad class of such problems. Here, we shall use some of these tools 
to study one particular problem of this class that arises in elasticity. 

In this paper we study a partially high contrast elastic composite whose 'inclu- 
sion' phase is disjoint and periodically distributed through the 'matrix' phase. The 
matrix phase is an arbitrary, generally heterogeneous material with a uniformly 
positive elasticity tensor and the inclusion phase is considered to be isotropic and 
'soft' in shear. Namely, the shear modulus for the inclusion material is chosen 
to be of the order e^, where e is the composite's periodicity size, while the bulk 
modulus remains uniformly positive. Such elastic inclusions can be called weakly 
compressible and elastic composites that are weakly compressible everywhere were 
studied in [16]. Composites with weakly compressible inclusions were first studied 
by Panasenko, via the method of asymptotic expansions, for e-independent body 
forces, c.f. [17, 18]. Therein, Panasenko showed the homogenised limit equations to 
be a system of equations depending only on the macroscopic spatial variable, i.e. 
not a two-scale system. A novelty of our work is that the body forces are allowed to 
be e-dependent, thus eventually allowing for a study of the corresponding spectral 
problem. The resulting homogenised limit equations are found to be genuinely two- 
scale, see Theorem 2.1, which reduce to the classical case for not only macroscopic 
body forces, as shown in [17, 18], but for a broader class of applied body forces, 
namely microscopically irrotational forces, see Corollary 2.2. We perform spectral 
analysis and show that the spectrum corresponding to the original problem con- 
verges in the sense of Hausdorff to the spectrum of the two-scale limit problem, see 
Theorem 2.3. Furthermore, due to the behaviour of the limit eigenfunctions pre- 
scribed by the weakly compressible condition, the two-scale limit spectral problem 
appears uncoupled in the macroscopic and microscopic problems. In the course of 
obtaining these results, we apply and develop further an appropriate modification 
of the two-scale convergence techniques to deal with the partial degeneracy, cf [15]. 

The layout of the paper is as follows: Section 2 contains the problem formulation 
and presentation of the main results. Section 3 reviews the necessary background 
material, in particular the method of two-scale convergence, its applications to 
homogenisation and further modifications to deal with partial degeneracies. Section 
4 and 5 are dedicated to the proofs of the main homogenisation theorem and the 
spectral convergence respectively. 
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2. Problem formulation and main results 




Q 




Figure 1. An example of a periodic composite body considered in this paper. 

We begin by introducing the notation used throughout this paper. U C R"^, d> 2, 
open and bounded, shall denote the domain occupied by the composite material, 
see Figure 1. Q = [0, 1)"^ is the periodic reference cell and shall consist of two 
disjoint regions: the 'inclusion' Q2, a subset of Q with smooth boundary T and the 
'matrix' Qi = Q\Q2- We assume that Qi is connected and Q2 does not intersect 
the boundary dQ, i.e. Q2 C (0, l)*^. We denote by F2 the Q-periodic extension of 
Q2 throughout R'^ and by eF2 the e contraction of F2, i.e. 

F2 := Ix : X = y + k iov some y G Q2 and some A; G Z > , 



and eF2 = {x : x/e G F2}. We denote by 01 and Of the inclusion phase and matrix 
phase respectively. That is, Q2 — ^^ ^^2 and fif = 0\ri2- For a Banach space X 
the spaces [X]'^ shall denote the space of vector- valued functions whose components 



belong to X, i.e. for u G [XY 



{ui,U2, . . ■ ,Ud} and mi,M2,- ■ ■ ,Ud G X. The 



spaces of matrix- valued functions [X] and fourth rank tensor- valued functions 
[X] shall be understood in a similar way. 

We shall consider the following problem: 



-V • ( C'e{u) )+au = f in f) 
u = on do,. 



(1) 



a > 0. The case a = is the elastostatic problem with body force /^ while a > 
is the resolvent problem, which is important for spectral analysis (cf. [11].) The 
underlying density function p(y) could be taken to be any positive matrix but is, 
for simplicity of deviation, assumed to be equal to unity. Here u G [^g (ri)] is the 

unknown displacement and e{u)ij = ^ ( ^ + ^ j is the infinitesimal strain tensor. 

Furthermore, C^(x) = C{x/e) for C{y) G [L'^iQ]]'^^'^^'^^'^, is the elasticity tensor 
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of the composite material which is considered to be Q-periodic, positive definite 
in the matrix and isotropic in the inclusion with Lame coefficients A ~ 0(1), 
H ~ 0{e'^). Explicitly C(y) is of the form: 

C(y) = C«(y)+e2cW(y), 
cijlqiv) = Xi{y)cljlg{y) + X2{y)Sijdpg, cljlg{y) = X2{y){Sip6jg + 6ig6jp). 

(We have set, for simplicity, A = 1, /i = ^.) Here 5 is the Kronecker delta symbol, Xi 
is the characteristic function oiQi, C^^^ G ^^^oof^Q^jdxdxdxd ^g ^g^j^gj^ ^q i,q symmetric 
and positive definite: 

^SlM = cf^M = ci%(y)^ Gg,(y)^P.%- > '^\v\', (3) 

for some v > 0, for all y € Qi, for all symmetric 77; /^(x) G [L^(il)]'^ is a prescribed 
externally applied 'body force'. 

We can see from (2)-(3) that C^^' is symmetric and non-negative while C^^'+C^^' 
is symmetric and positive definite, i.e. 



^S.(y) = ^]L(y) = ^S.(2^)' for Z = 0,1. (4) 

cijUy)vp,mj > 0, {cijigiy) + C^^lgiy)) ^p,% > Hr/|^ 

for some z^ > 0, for all y £ Q, for all symmetric rj. In (4) and henceforth summation 
is implied with respect to repeated indices. We now state our first main result. 

Theorem 2.1 : Let /^(x) weakly (strongly) two-scale converge to f{x,y) as e ^ 
0. Then the sequence u^ , of solutions to (1), weakly (strongly) two-scale converges 
to u^ix,y) = u{x)+v{x,y) as e ^ 0, where {u,v) £ [H^{n)f x [L^{n; H^{Q2))f 
is the unique solution to 

-V • (c^^™"Vm(x)) + au{x) + a{v){x) = {f){x) m ft, (5) 

-Ayv{x, y) + av{x, y) + au{x) = f{x, y) + Vypix, y) in Q2 

Vy ■ v{x, y) = in Q2 (6) 

v{x,y) = on dQ2, 

where p G Lp' {ft; H^ {Q2)) is unknown. (•) denotes the mean value over Q, i.e. 
(/)(x) = fQf{x,y) dy. C ™ is the constant coefficient positive definite tensor 
given by 



^'°"-/ ^^.'V?^) (m.. + ^J'' 



C^ = ClLiy) 6prS,s + ^^ dy. (7) 



Here Nrs = {N^^, ^rs^ • • • ) ^rs) "^^ ^ Q-periodic solution to the degenerate cell prob- 
lem 

-Vy(c^^\y){er'^es + VyNrs{y)))=0 inQ. (8) 

The solution Nrs to the degenerate cell problem (8) is not unique as non-trivial 
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solutions to the homogeneous problem exist due to the degeneracy of C^^' . Such 
homogeneous solutions form a non-trivial subspace of [H],{Q)] and are charac- 
terised by satisfying the condition C^^'{y)'VyV = m Q. However, we can see 
from (7) that (7^^°™ remains unchanged if we add any homogeneous solution of the 
degenerate cell problem to N; from this we can conclude that even though the 
solutions to the degenerate cell problem are not unique, (7'^°™ is unique. Further- 
more (7^^°™ is positive definite, which can be seen by noticing from the variational 
form characterisation of the tensor (7^^°™ that it is "more positive" than the effec- 
tive tensor corresponding to the classical homogenisation of a perforated domain 
problem, which is indeed positive, see e.g. [19]. 

Let us now consider problem (6). One notices that for microscopically irrotational 
body forces, namely forces of the form f{x, y) = fo{x) + Vyfi{x, y), one can choose 
p in (6) in such a way as to 'absorb' the forcing term, i.e. by setting p{x, y) = —y ■ 
(/o(x) — au{x)) — fi{x, y) +p{x, y). As a result, (6) is reduced to the homogeneous 
Stokes problem: find v,p, such that 

-Ayv{x,y) + av{x,y) = Vyp{x,y) in Q2 
\7y ■ v{x,y) =0 in Q2 
v{x,y) = on dQ2, 

which is well known to have only the trivial solution v = 0, see e.g. [20]. This 
shows that for microscopically irrotational body forces the limit problem (9) is 
independent of the microscopic variable y, i.e. the limit solution u^{x,y) = u{x) 
has no microscopic oscillations. Furthermore, in this case u^ ^ u strongly in L^ 
when f^{x) -^ f{x, y). These striking results are very different to the double poros- 
ity case where that the homogenised limit is of a genuine two-scale nature, for a 
general body force /*", with the limit function depending on both the macroscopic 
and microscopic variable. The reason for the limit function having no microscopic 
oscillations, for the described body force, is precisely due to the form of the partial 
degeneracy: the inclusion phase is isotropic with Lame coefficient A ~ 0(1), which 
means that in the asymptotic limit e — )■ the inclusion phase is microscopically in- 
compressible, i.e. Vj^ -v = 0, and, as a consequence, for microscopically irrotational 
body forces no deformations will occur on the microscopic scale. The following 
corollary states these observations in a precise form. 

Corollary 2.2: Let f^{x) strongly two-scale converge to f(x,y) = /o(a^) + 
Vy/i(x,y) for given, sufficiently regular, /c/i- Then the sequence u^ , of solu- 
tions to (1), strongly converges to u{x) in L^(r2) as e — )■ 0, where u G [H^Y ^-^ ^^^ 
unique solution to 



V 



(c^°'^Vn(x)) + au{x) = {f){x) in Vi. (9) 



Postponing a precise definition until Section 5, let A^ and A^ be the self-adjoint 
operators corresponding to problems (1) and (5)-(6) respectively. We shall now 
state the second main result of this paper. 

Theorem 2.3 : The spectrum of A^ , ct{A^), converges in the sense of Hausdorff 
to the spectrum of A^, a{A^). That is 

(1) For every A" € f^{A^) there exists A^ G o'l^^) such that A^ — t- A''. 

(2) If there exists A^ G a{A^) such that A= -^ A°, then A° G criA^)- 

Theorem 2.3, whose proof is given in Section 5, tells us that if we wish to study 
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the limit behaviour of the eigenvalues of A*" as e — )■ then it is sufficient to study the 
spectrum of the limit problem a{A^). To this end consider (7{A^) and let {X^,u^) 
be an eigenvalue-eigenfunction pair of A^. Then u^(x,y) = u{x) + v{x,y) satisfies 

-V • (c^'°'"Vn(x)) = A%(x) + A°(t;)(a;) in ft, (10) 

-Ayv{x,y) = X'^v{x,y) + Vyp{x,y) m Q2 
Vyv{x,y) = iuQs (11) 

v{x, y) = on dQ2- 

Note that the equations (lO)-(ll) are not coupled: their uncoupled nature is due 
to the fact that X^u{x) originally present on the right hand side of (11) can be 
absorbed by p{x,y), i.e. u{x) = Vyq{x,y) for q{x,y) = y • u{x,y) for y G Q2. The 
main consequence of this uncoupling is that the spectrum of A^ is simply the union 
of the spectra corresponding to problem (10) and (11) respectively, i.e. 

Corollary 2.4: The spectrum of the homogenised lim.it operator A^ , a{A^), has 
the following representation: 

cr(^°) = {A„ I n G N} U {;[/„ \ m G N}, 

where A„ satisfy, for some non-trivial Un G [i:fQ(r2)] , 

-V • (c^°'^VUn{x)^ = XnUnix) in Q, 

and fim satisfy, for some non-trivial Vm G [Hq{Q2)]'^, Pm G H^{Q2), 

-AyVmiy) = fJ"mVm{y) + "^ yPmiv) in Q2, 

Vy • Vmiy) =0 in Q2. 

Remark 1 : In the case of il = R the operator Lu := —V • (C ™Vu(a;)) has 
absolutely continuous spectrum coinciding with the positive real line. The spectrum 
of A^ is still the union of the spectra for the operators defined by (10) and (11), 
which implies that the spectrum of A^ contains no gaps; o'{A^) = [0, oo) and 
contains eigenvalues of infinite multiplicity /U„ corresponding to the eigenvalues of 
the Stokes problem on the inclusion Q2. Furthermore, the spectral convergence 
result. Theorem 2.3, can be shown to hold in this case too. 

3. On two-scale convergence and its modifications for partial degeneracies 

Let us review the concept of two-scale convergence and some of its properties that 
are useful in the homogenisation of second order PDEs. For a full account of two- 
scale convergence and its application to homogenisation theory see e.g. [8~11]. 

3.1. Two-scale convergence and strong two-scale resolvent convergence 

Definition 3.1: (Weak two-scale convergence.) Let u^ be a bounded se- 
quence in L^(r2). We say u^ (weakly) two-scale converges to u^ G L'^{Q, x Q), 
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denoted by u^ ^ u°, if for all (/. G C^{n- C^{Q)), 



u^(x)0 (x, I) dx — > u (x, y)0(x,y) dxdy 

n Jn Jq 



as e — )■ 0. 



The following properties of two-scale convergence are of particular importance 
to us: 

Lemma 3.2 Properties of (w^eak) t^vo-scale convergence 

(i) If u^ is bounded in L'^{fl) then there exists u" G L^(r2 x Q) and a subsequence 

n^' such that u^' -^ u^ . 
(a) If u^ ^ u^ then n^ converges to JqU^ dy weakly in L^(fi). 

(Hi) If u^ -^ u^ and a{y) € L'^{Q) then a(x/e)u^{x) -^ a{y)u^{x,y). 
The following lemma is of particular importance to high contrast homogenisation: 

Lemma 3.3: Let u'^ G H^{^) such that 

IP Wv^isi) - '-'' e\\\/u ||2,2(f^) ^ o, 

for some constant C independent of e. Then, there exist u^{x,y) G L? {yi; H\{Q)) 
such that, up to extracting a subsequence in e which we do not relabel, 

u'^u\x,y) (12) 

eVu^ ^VyU^{x,y). (13) 

In this paper we shall consider a sequence of non-negative self-adjoint operators 
A^ and wish to study the limit behaviour of their resolvents, i.e. we wish to study 
the 'resolvent' problem: for fixed a > 0, 

A^u + au = f, 

as e — )■ 0. In homogenisation theory we often find that u^ := (A^ + a)^^f G L'^{Q) 
two-scale converges to some vP{x,y) G L'^{Q x Q) where 

A\^ + an° = /, 

for a self-adjoint operator A^. The notion of strong resolvent convergence is not 
applicable here since the limiting operator [A^ + a)^^ is defined on a subspace of 
L^(J7 X Q) while {A"^ + a)~^ is defined on L^(r2). We instead use the notion of 
strong two-scale resolvent convergence which we shall now recap, see [10, 11] for 
more details. 

Definition 3.4: (Strong two-scale convergence.) A sequence u'^ G L^($7) is 
said to strongly two-scale converge to u ^ L^{Q, x Q), denoted n*" -^ u, ii u^ ^ u 
and 



u^(x)f^(x) dx — )■ / / u{x,y)v{x,y) dxdy, for all v^ -^ v. 
n Jn Jq 
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Equivalently: u^ ^ u if, and only if, 



u^ ^ u, &: lim / (n^) (x) dx = u (x, y) dxdy. 

^^0 Jn Jn Jq 

Definition 3.5: Let A^, A^ be non-negative self-adjoint operators on L^($7) and 
a closed linear subspace H of L^($7 x Q) respectively. Then we say A^ strongly 
two-scale resolvent converges to A^, denoted A^ A A^ if for some a > 

{A^ + a)-^f A (A° + ay^Pf whenever f ^ f eL^{nxQ), 

where P : L^{Q x Q) ^ H is the orthogonal projection onto H. 

Remark 1 : Since the resolvent set is an open subset of C and the resolvent is 
an analytic operator- valued function on C it is sufficient to test, as in the case of 
strong resolvent convergence, strong two-scale resolvent convergence for a single 
a > in the resolvent set, say a = 1. 

If A^ strongly two-scale resolvent converges to A^ then the limit spectrum cr{A^) 
is always contained in the limiting spectrum linig^o o'(A^). That is 

Lemma 3.6: Let Al" , J^ be non-negative self-adjoint operators on L^(r2) and 
the subspace H C L'^{0, x Q) respectively, let A^ ^ A'^ . Then for all A° G o-(^°) 
there exists X^ £ cr(^^) such that A^ — )■ A*^ as e ^ 0. 

In general, it is not the case that the reverse inclusion holds. However if true, i.e. 
if A^ 6 cr(^^) such that A"^ — ^ A° as e — ^ imphes A° G a[A^), the proof is more 
difficult and problem specific. Such proofs usually require establishing, by separate 
means, a version of 'two-scale spectral compactness', see Section 5. 



3.2. Modification of two-scale convergence for partial degeneracies 

An important auxiliary problem in the homogenisation of problem (1) is the so 
called "degenerate cell problem": find v G [-ffi(Q)]'^ such that 

-Vy(cW{y)Vyv)=F, (14) 



where F G Hn (Q) is given. Hj, (Q) is taken to be the dual space of [H],{Q)]'^ and 
we denote, for F G H^ (Q), by {F,w) the duality action of F on tt; G [-ffl(Q)]'^- 

The weak formulation of (14) is: find v G H],{Q) such that 

C^'Hy)Vyv{y)-Vyw{y)dy ={F,w), V u; G [i/^(Q)]' . (15) 

Introducing the space 

V:={u€ [HliQ)f : C^'Hy)VyU = o} , (16) 

we see that if v solves (14), equivalently (15), then F must necessarily satisfy 
{F, u) = for all u £ V. This turns out, and was first shown by I.V. Kamotski 
and V.P. Smyshlyaev in [15], to be a sufficient criterion for solvability of (14) 
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when making the following Key assumption on the degeneracy: There exists 

d 



a constant C > such that for all v G 



H^iQ) 



there exists vi € V with 



,rl < C 



c«(y)v,^ 



The condition (17) can be equivalently re- written as 



(17) 



l-fV^l'llrrrl /^^l'' < C 



[Hl{Q)Y 



c^'\y)Vyv 



(18) 



where Py^ is the orthogonal projector in H],{Q) on V , the orthogonal com- 
plement to V. The assumption (17) does not depend on the choice of an equivalent 
norm in H],. The assumption (17) holds for most of the previously studied cases, 
and has to be checked by separate means, see Section 4.2. 

Lemma 3.7: 

r T^ 

(i) Under the assumption (17), the problem (15) is solvable in H],(Q) if and 
only if 



{F,w) = 0, yw ev. 



(19) 



When (19) does hold, the problem (15) is uniquely solvable in V . 

(a) For any solution v and any vi £ V , v + vi is also a solution. Conversely, any 

two solutions can only differ by a vi £ V. 

For completeness we will present the proof to Lemma 3.7 which was first proved 
in a greater generality by I.V. Kamotski and V.P. Smyshlyaev, cf [15]. 

Proof: (i) Let f be a solution of (15) and let w £ V. Then, using the symmetry 
ofC(i) and (16), 

{F,w) = [ C^'\y)Vyv{y)-Vyw{y)dy = [ Vyv{y) ■ d^\y)Vyw{y) dy = 
JQ Jq 

(20) 

yielding (19). Conversely let (19) hold, and seek v G H\{Q) solving (15). By 

(20), the identity (15) is automatically held for all w in V , therefore it is sufficient 
to verify it for all w G V-^ . Seek v also in V-^ . Show that then, in the Hilbert space 

r -id 

H := V with the inherited H\{Q) norm || • \\h, the problem (15) satisfies 



the conditions of the Lax-Milgram lemma, see for example [21]. Namely, first the 
bilinear form 



B[vM := / C^^\y)Vyv{y)-Vyw{y)dy 



is shown to be bounded in H, i.e. with some C > 0, 



i?[f , w\ 



< C \\v\\h\\w\\h, yv, w G h. 



This follows from (2) where C(2) g [L°°((5)]'^^^^'^^^. We will now show that the 
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10 

form B is coercive, i.e. for some u > 0, 

B[v,v] > iy\\v\\jj, yv £ V^. 

We have 



B[v,v]:= d^\y)Vyv{y)-Vyv{y)dy= (cW(y)) V 



1/2 



> 



c 



c^^\y)'^yv 



> U \\v\ 



H- 



1 /9 

In the last two inequahties we have used, the boundedness of [C^^'^ and (18). 

Therefore, by the Lax-Milgram lemma, there exists a unique solution to the 
problem 

veV^: B[v,w] ={F,w), 'iw e V^ , 

and hence to (15). 

(ii) If V solves (15) and vi ^V then C^^> {y)V yVi{y) = and hence v + vi also 
solves (15). Assuming further v^"^' and v^"^' both solve (15), v := v^^' — v^^> solves 
(15) with F = 0, and then set w = v. As a, result. 



C^^\y)Vyv{y)-Vyv{y)dy 



'Q 

(1)aV2 



1/2. 



c«(y))^ V 



implying (^C^^'^ Vj^f = and hence C^^'VyV = 0, i.e. v ^ V. 



D 



4. Homogenisation 

4.1. Proof of Theorem 2.1 

The weak form of (1) is stated as follows: Find u'^ € [i7g(r2)] such that 



C(i) (f ) e{u') ■ e{(t)) + e^cW (|) e(n^) • e(0) + au^ ■ (j) = / f ■ (f) 



^{Hlm''. (21) 



For any fixed e > 0, the associated quadratic form 



A'{u, v) := / C(i) (f ) e{u) ■ e(t;) + e^cW (|) e(u) • e{v) + au-v, (22) 
is coercive, i.e. there exists C > such that 

This is established by noting C'^^' + C^^' is positive definite, see (4), and an appli- 
cation of the standard Korn's inequality, see e.g [19]. By the Lax-Milgram lemma. 
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c.f. [21], coercivity ensures the existence and uniqueness of n^ G [Hq{Q)] for fixed 
a > 0. We wish to study how u'^ behaves as e tends to zero. Substituting (p = u^ 
into (21) and using (4) we see that there exists a constant C independent of e such 
that 



\L^{n) 



<C\ 



\L^{n) ' 



|eV-u'^||^2(f^) < C||/'^||^2(f^) 



(C(i)(x/e))i/2Vn^ 



L^{n) 



<C\ 



\L^{n) 



(23) 
(24) 

(25) 



Note that for a = we require the vaHdity of the foUowing Poincare type inequahty: 
There exists a constant C > independent of e such that for ah u G Hq{Q) 



l^llL2(n) ^ ^ 



C(i)(x/e)Vn 



LHni) 



+ e^ 



C^^\x/e)Vu 



LH^i) 



which indeed holds, the proof of which is given in Section 5.1, Lemma 5.3. The 
uniform bounds (23)-(25), and the two-scale compactness result, see Lemma 3.2 
(i), imply that, for uniformly bounded f^, these sequences have two-scale conver- 
gent subsequences and the behaviour of their two-scale limits is described in the 
following theorem: 

Theorem 4.1 : There exists u{x) G [H^{n)f, v{x,y) G L'^i^; [^^(^2)]'^), Vj^ • 
v{x,y) = 0, such that, up to a subsequence in e (which we do not relabel), 

u^ -^ u{x) + v{x, y), 
eVu'^ ^ Vyv{x,y), 
{C^^Y\x/e)Vu' ^ (C«)'/'(y) [VX^) + VX(x,y)] , 

where u^{x,y) G [L"^ {Q; H},{Q))]'^ is a solution to 



Vy ■ {^C^'Hy)VyuHx,y)) = V, • (^C^'\y)V M^)) ■ 



(26) 



Proof: 

Step 1: By (23)-(24) and Lemma 3.3, up to a subsequence in e, which we do not 

relabel, u^ and eVu^ two-scale converge to u^ and VyU^ respectively. 

Let us show u^ G [L'^{Q;V)]'^, for V given by (16). Since eVu^ -^ ^yU^, by 
Definition 3.1 and Lemma 3.2 (iii). 



\l/2 



(C«(f)) Vn-(x).0(x,f) 



dx 



1/2 



eVu^{x)- ("cW [ff\ (t){x,l) dx 



n JQ 



Vj,nO(x,y)- C«(y) cP{x,y) dydx 



1/2 



n JQ 



1/2 



C(i) (y) Vyu^{x,y) • ^{x,y) dydx, \/^ G [C^ {Q; C^ {Q)f ■ 
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By (25), {C^'^\x/e))^/'^Vu''{x) is bounded in L'^{n) so 

I e(cW (f) )^^Vu=(x) • (j) (x, f) dx ^ 0. 
Therefore 

/ / {d^\y)f'^Vyu\x,y) ■ <Pix,y) dydx = 0. 
Jn JQ 

This impHes, since the functions (j){x,y) are dense in [Lp'{Q. x Q)] , that 

{C^^\y))'^/^VyU^{x, y) = a.e. x G f^. (27) 

Premultiplying (27) by {C^'^\y)Y/'^ shows vP{x,y) G ^^(ri;^). For a.e. x e Q., 
uP{x,-) G V, and by (50) below, we see u^{x,y) = u{x) + v{x,y) for some u G 
[L\n)]^, V £ [L^n;H^{Q2))]^ with Vj, • t; = 0. 

Inequahty (25) and Lemma 3.2 (i) imply that there exists ^^ G [L^(J1 x 
Qjjdxd such that, up to a subsequence in e which we do not relabel, 
{C^^'{x/e)y''^'Vu'^{x) -^ ^^{x,y). Introducing the space 



W^ := I* G [L2(Q)]'^x'^ : V, • (^(c«(y))'^' ^(y)^ = in i/#HQ)} 



(28) 



we will show that ^" G L^{n;W). Take in (21) 0(x) = (^=(x) = e$(x, f) for 

any <I>(x, y) G C^ i il; C^{Q) I . Passing then to the limit in (21) we notice, via 

(23)-(25), that the limit of each term but the first one on the left hand-side of (21) 
is zero, and therefore 

lim / C(i) (-^ Vn^(x) • eV$ (x, -) dx 

1 /9 

C(i)(y)) e{x,y) ■ VyHx,y)dydx,= 0. 




inJQ 
The density of '^(x, y) implies then that, for a.e. x, 

Vy • ( (c«(y))'^'e°(x,y)) = in H-\Q). 

This yields £,^{x,y) G Lp' (il; W), see (28). Hence, we have shown 

u^ -^ u{x) + v{x,y), 
eVu^ ^ Vyv{x,y), 
(CW)i/2(x/e)Vn^^C0(x,y), 

for some u G [L'^i^)]'^, v{x,y) G ^^(Jl; [//^(Qa)]'^), V^^ • f = and C^{x,y) G 

L2(17;VF). 
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Step 2: Let us now show that u{x) G [Hq{Q)]'^. For ^ £ C°°{0,; W) we see via 
integration by parts 

/ (C^^^ (f) )^^^Vn^(a^) • ^ (^> f) dx = - f u' ■ (v^ ■ (c^^'> (f) V^^^ {x, f) 

and passing to the two-scale Hmit indicates u^{x, y) and 6.^{x, y) are related by the 
following expression 

e{x,y)-^{x,y)Aydx = - f f u\x,y) -V ^ ■ ({C^^\y)Y'^^{x,y)) dydx 
qJq JqJq ^ ' 

(29) 

Now, for fixed ip G [(7°^(ri)] we take in (29) test functions ^(^ given by Lemma 
4.2 at the end of this section. Then, (29) states, via (42), 

l{ip) = - u{x) -V^-if dx, (30) 

Jo. 

where 

^{'f)-= I / ^o(x,t/) -^^(x,?/) dydx, 

Jo. JQ 

defines an L^ bounded linear form on (^ G [C°^(i7)]'^^'^, see (43). Therefore, by the 
Riesz representation theorem, (30) implies 

i{4>) = V ■ If dx, 
Jn 

for some v G [L^{n)f'''^ and hence, via (30), u G [H^{n)f. 

Step 3: It remains to show ^^{x, y) = {C^'^'Y''^{y)^xu{x) + \IyU^{x,y)] for some 
v} given by (26). For a.e. x G $7, let ^-'^(x, •) G [H\{Q)Y be a solution of (26). Note 
that, if the key assumption (17) holds, such a solution exists by Lemma 3.7 since, 
forF:=Vj,-(C(i)(y)V,n(x)), 



(F, v) = - C(i)(y)V,n(x) .Vyvdy = - V,n(x) • C^'Hy)VyV dy = 
JQ JQ 

VwG V. 

The key assumption (17) does indeed hold and the proof of this fact is the subject 
of Section 4.2. Setting ^(x, y) := (C(^))^/2[Vj,n(x) + VyU^{x, y)], we note C is well- 
defined, since n^(x,-) is unique up to a function in V. Furthermore, ^(x,y) G 
[L'^{n;W)f, from (26), with 

ax,y)-^{x,y)dydx = - [ f u{x) -V^- ({C'^^\y)f'^^{x,y)) dydx 
nJQ JnJQ^ ' 

V^G C°^((7;VF). 
(31) 
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The latter directly follows via integration by parts and (28). Now let us show that 
£,oix,y) = ^{x,y) for a.e. a; G $7: since both ^o(a^5 •) and S,{x, •) belong to W, it is 
sufficient to show S,o{x, •) — ^(x, •) _L W. To this end, (29) and (31) imply that 




{^o{x, y) - C{x, y)) ■ ^(x, y) dydx 
n JQ 



n JQ 



z;(x,y)-V,-((C«(y))^/2vI/(x,y))d|/dx V^ G C~([7; VF). (32) 




Now the result follows if the right hand side of (32) is zero. Since C^{Q.;V) is 
dense in L'^{Q.; V) it is sufficient to show that for a fixed 4> £ [C^i^; H^iQ2)riV)f 

(C(i)(y))^/2V^0(x, y) • ^ dydx = V^ G C°^(J1; W) (33) 

The latter can be seen by considering a function A^ G [H},{Q)]'^ such that N{y) = y 
for y G Q2- An example of such a function would be an appropriate periodic 
extension of f{y) = X2{y)y, which is possible since Q2 C (0, l)'^. Now for given (p 
define an auxiliary function ^{x,y) := 'Vx4>ix,y) ■ N{y). For a.e. x £ 0,, $(x,-) G 
[H^iQ]]'^, and for any * G VF 

= (-V, • ((C(i) (y))'/^^) , Hx, y))= I ^{y) ■ (C^^) iy))'/^VyMx, y) dy 

*(y) • {C^^\y)f/^Vy{y.4>{x,y) ■ N{y)) dy 

vl/(y) • (C«(y))V2v,(V,0(x,y) -y) dy 

'2 

vI/(y)-(C«(y))V2v,</,(x,y)dy, 



where the last equality is due to the fact that, for (/)(x, •) G F, in components, 

(C«(y))i/2v,(V.0(x,y) -y))^^. = {C^'\y))]il{h,..ys),y, 

= (.C^^Hy))]jpqi(I^P,x.y,ys + (l>p,xjsq) 

= ic^'Hy))]Sqi^P,x.y,ys) HC^'Hy))]Sqi^p,x:. 

^ ^ ' 

=0 

(C«(y))i/2v,(A(x,y) 

Here, we have used the notation , z to denote differentiation with respect to the 
variable z, for example fpp^x, '■= t^- D 

Using Theorem 4.1 we are now able to pass to the limit as e — )■ in (21) and 
prove Theorem 2.1. 

Proof: of Theorem 2.1. 
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,0^ 



Step 1: For fixed (t){x,y) G Cq''{Q.;V), employing tlie test function (jr{x) : 
(x, |) in (21) gives 

/ C(i) (f ) V^^ • V,A (x, f ) + eC(0) (f ) Vu' ■ [eV,0 (x, f ) + Vy<t> (x, f )] + 
Jo. 



+ / a^^-0(x,f)= / r-0(x,f). (34) 
Jo. Jn 



Passing in (34) to the limit e — )■ 0, using Theorem 4.1, gives 
C(i)(y) [VMx)+^yuHx,y)] • V^0(x,y) dydx+ 




n JQ 

+ 11 C^^Hy)^yvix,y) ■Vy(l){x,y) + a{u{x) + v{x,y)) ■ (l){x,y) dydx 
Jn JQ 

= f f f{x,y)<l)ix,y)dydx, y<l,{x,y) € [C^in;V)f. (35) 
Jn JQ 

Step 2: Choosing in (35) (p{x,y) = (p{x), ip G Cq°{Q,), gives 

C^^\y){Vxu{x) +VyU^{x,y)) -Vx^ix) + a{u{x) + v{x,y)) ip{x) dydx 
n JQ 

= fix,y)(p{x) dydx, 

JflJQ 

which is the variational form for 
-V ■ ((c^^\VxU + VyU^)) (x)) +a{u{x)+ < V > (x)) =< f > (x) in Jl. (36) 




Setting Up = Nrs{y)-^{x) and substituting into (26) and (36) gives equations (5) 
and (7)- (8). Let us show that (7*^°™ is strictly positive. For any symmetric rj G R'^^'^ 
we can show, as in the case of the perforated domain (see e.g. [22]), 0^°"^!^ ■ rj has 
the following variational representation 

C^^^^^ . r/ = inf / C7(i) {y) {t] + Vyw) [t] + V^,^) dy. 

we[c^{Q)YjQ 

By (2) it is easy to see that 

C^°'^r/-r/> inf / C^^\y) (r] + V yw) (f] + V yw) dy 

■'^(^[c$(Q)YJq^ 

= Cr]-r], 

where C is the homogenised tensor for the perforated linear elasticity problem 
which is well known to be strictly positive, see e.g. [22]. Therefore (7^^°™ is strictly 
positive. 

Step 3: Now choose in (35) the test functions of the form 0o(x,|/) = ^(x)0(y), 
^ G C(f (f]), (j) G C^{Q2) with V • (/) = 0. Due to (33) and the fact Co{x,y) = 
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(CW(2/))^/^ [V^u(x) + VX(2;,y)] G [L'^{^;W)Y, we have 

/ / [C^^Kv)) [^xu{x) + Vyu\x, y)] V^ij{x)4>{y) dydx = 
Jn JQ 

Furthermore 



C''^\y)'^yv{x,y)-Vy(l){y) dy(lx= I / {5ip5jq + 5iq5jp)vp^y^<j)i^y^ dydx 
nJQ2 JnJQ2 




{vp,y,(t>p,y, + Vp,y^<pq,y^) dydx = / / (i^p.y, (/>p,j/, + Vp,y^<f>q,yJ dydx 
nJQ2 Jn. " 



VyV • VyCp dydx. 
inJQ2 



Therefore, 

Vyv{x, y) ■ Vy(j){y)il){x) + a {u{x) + v{x, y)) ■ il){x)4>{y) dydx 




nJQ2 

= 11 f{x,y)'4){x)(t){y) dydx, 

JnJQ2 

which is the variational form for (6). 

Step 4: It remains to show that ii f^ ^ f then u'^ A n°. This proof is similar to 
the double porosity case, see e.g. [10], and we will present it here for completeness. 
Let z^ £ [Ho{n)f be the solution to 

-V-{C^Vz^) + az^ = u^. (37) 

Then, by the above arguments, z^ -^ z^ where z^{x, y) = z{x)+w{x, y) is a solution 
to 



wll(y) [^w^qs + N?s,q{y)] Zr,s{x)(t)i,x^{x,y) dydx+ 
Q..JQ 

+ 11 C^°Hy)^ywix,y) ■Vy(j){x,y) + a{z{x) + w{x,y)) ■ (f){x,y) dydx 
Jn JQ 

= [ [ u\x,y)<P{x,y) dydx, y<Pix,y) G [C^in;V)f. (38) 
Jn JQ 

Setting (j) = z^ in (35) and (p = u^ in (38) shows 

u^ix,y)\ dydx = / / f{x,y)z^{x,y) dydx. (39) 

nJQ JQ.JQ 

Similarly, the variational forms for (37) and (1) show 

M=(x)|2 dx= I fix) ■ z'{x) dx. (40) 

n Jn 
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Using the assumption f^ ^ f^, passing to the hmit e — )■ in (40), via (39), gives 

hm / \u^\ dx = hm / f^{x) • z^{x) dx = / f{x,y) ■ z {x,y) dydx 
£^0 J^ e-^o J^ J^ Jq 

\u {x,y)\ dydx. 



n JQ 

Hence, by Definition 3.4, u^ A u^. 

Finally, observe that the uniqueness to the solution of (5)-(6) follows in a stan- 
dard way by setting in (35) (j) = u^ and / = 0. D 

Proof: of Corollary 2.2. The equation for the microscopic deformations v{x,y) 
is given by 




Vyv{x, y) ■ Vy'ip{y)(p{x) + a ( u{x) + v{x, y) ) • (/>(x)V'(y) dydx 
nJQ2 



f{x,y)4>{x)'ip{y) dydx, (41) 
nJQ2 

for all (J) G C'(f (r^), V G C^iQ2) with V • V = 0. For an externally apphed body 
ioice f%x) = f{x,x/e), where /(x,y) = foix)+Vyfi{x,y), such that f A f{x,y) 
as e — )• we find, by taking into account that for any constant vector field c we 
can have the representation c = Vy{c • y) in L'^{Q2), that 



/(x, y) ■ 4>{x)ij{y) dydx = / / [foix] + Vj//i(x, y)] (t>{x)ip{y) dydx 
n JQ2 Jn 



i^y (M^) ■ y) + Vj;/i(x, y)] (t){x)^{y) dydx 
nJQ2 



[fo{x) ■ y + fi{x,y)] 4>{x)\7y ■ ip{y) dydx = 0, 
IQJQ2 

for all ^p £ C^{Q2) with V •-(/; = 0. Similarly, J^Jq u{x) ■ (l){x)ip{y) dydx = 0. 
Therefore equation (41) becomes 



Vyv{x, y) ■ \7yip{y)(j){x) + av{x, y) ■ (t){x)ip{y) dydx = 0. 
'nJQ2 

This implies that for a.e. x in fi, v{x, •) is a weak solution of the homogeneous 
Stokes problem which is well known to have only the trivial solution ^(x, y) = 0, 
see e.g. [20]. Furthermore, as u^ A u{x), u^ ^ u strongly in L^. D 

Let us end the section with the proof of the Lemma used to show the regularity 
of the limit function n(x). 

Lemma 4.2: For fixed if G [C'^{VL)Y, there exists ^{x,y) £ (7°°(17;Ty), such 
that ^ = in Q2, and 

V.-(/. = V,-('y (C«(y))'^'*(x,y)dy') (42) 
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Furthermore, there exists a constant C independent of (p such that 

1/2 / /■ /■ \ 1/2 

('" 

Proof: Let us introduce the tensor 

dNL 



I h-^ <c(f \ipA ( f f lei' dydx) , Ve G L\n; W). (43) 
Jn JQ \Jn / \Jn Jo / 



a 



tjrs •- 



C^I,(JJ) { ^W^s + 



for Nrs G [Hi{Qi)]'^ being the solutions to the ceU problem 

-Vy(c^^^VyNrs)=0, mQi, 



(2) dN? 



-C-frs^j, on dQ2. 



(44) 



The tensor C is well known to be the positive homogenised elasticity tensor for a 
perforated elastic body with elasticity tensor C^'^\ see e.g. [22]. 

For a fixed if £ [C°^(rj)]'^, let ip°{x) G [H^{^)]'^ be the unique solution to 



V^ • ((7V^(^°) = V^ • V3. 



(45) 



Such a (p^ exists by the positivity of C. Furthermore, ip'^ G [C°^(r2)]'^, (p^\dn = 0. 
Next, take ip^ G C°^(ri; H^iQi)) to be a solution to 



-Vy ■ (C(2) [Vy^\x, y) + V,/(x)] j = 0, in Qi, 
(7(2)Vj,(^^(x, y)-n = -C^'^^V^p^{x) ■ n, on ^Qs- 



(46) 



Setting 



1/2 



^(x, y) := xi{y) [C^^\y)) \^y^\x, y) + V,/(x) 



(47) 



by construction, \I' G C°^(r2;VF), see (28), ^{x,y) = for y G Qi- Furthermore, 
p>l{x,y) = N^s:^, for Nrs solving (44), and 



Qi 



1/2 



cW(y) ^{x,y)dy 



ClLiy) ^rS,s + 



cm 

dVq 



dy 






(48) 

Therefore, (45) and (48) imply (42). It remains to prove (43). For fixed ^ G 
L'^{^]W) 




^ • ^ dydx = 
IQJQ JnJQ 



1/2, 



e-(c(2)(y)) V.p'dydx 



< I I \VxV" 



1/2 



f f C'^^k • ^ dydx 
Jn Jqi 



1/2 
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Now inequality (43) follows from (3) and observing that 



<.(|jV./|^y'7/'|v.|^ 



1/2 / /■ \ 1/2 

n 



where the first inequality is due to the positivity of C and then recalling (45). D 

4.2. Proof of main condition (17) 

In this section we prove the validity of the key assumption (17) for the tensor C^^> 
given by (2)-(3). That is we wish to show: there exists a constant c > such that 
for any uG [H\{Q)Y 

\\PvM\]i^{Q) < c (l|e(u)||i2(Q^) + ||V • u\\l,^Q^)^ . (49) 

Here Py± is the projection on to V-^, the orthogonal complement of 

V = {ve [HliQ)]'' : v{y) = k + X2{y)w{y) for some fc G R^ u; G [^^0(^2)]^ 

VyW = Q}, (50) 

see (16). The inequality (49) holds trivially for any element of V and need only be 
validated for elements of V . Moreover, it is sufficient to show inequality (49) for 
the following equivalent ifi((5)-norm 



\H ■- 



u dy 
Qi 



2 

2 



+ / |Vnr dy, (51) 



using the convention f u dy = (^f ui dy, . . . , f u^ dy) for u = (ui, . . . ,Ud)- The 
norm (51) is induced by the following inner product 

{u,v)h-=( udy]-( V dy] + Vu-Vvdy, (52) 

V^Qi / \Jqi / Jq 

and by definition, w G V if 

(/ w dy] ( vdy]+ VwVvdy = Vv G F. (53) 

As constant vectors are in V, see (50), Jq w dy = and (49) easily follows from 
(51) and the following result: 

Lemma 4.3: There exists a constant c > such that 

\Vw\^dy<c( \e{w)\^dy+ |V • wp dy ) Vu; G F"^. (54) 

Jq \Jqi Jq2 ) 
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Proof: Functions (f) G [C^{Q2yf with V • = belong to V, and by (53) w ^V^ 
if 



{)={w,(t))H= / Vtt;- V(/) dy = -(Au;,0), (55) 

where (Aw, cj)) denotes the action of the distribution Aw on the test function cj). 
Equation (55) states the distribution Aw is orthogonal to all divergent free test 
functions in Q2. It is well known, see e.g. [20, Proposition 1.1, p. 14], that such 
distributions are potentials, i.e. Aw = VV' for some distribution ij: and, since 
w G [H^{Q2)Y, ip G L'^{Q2)- Therefore, we see that w G V-^ if, and only if. 

Aw = VV' in Q2 for some V' G L^{Q2) , / w dy = 0. 

For fixed w G V , let w be its harmonic extension, see Lemma 4.5 below. Denote 
wi := w — w. Evidently, w = w + wi, wi G [Hq{Q2)]'^, Jq w dy = 0, and 



w = w in Qi, Aw = in (52, 
wi = in Qi, Awi = Vif in Q2, 



where f G L'^{Q2)- Since 



|Vw| dy < / |Vw| dy + |Vwi| dy. 



and 

/ |V-wi|^</ |V-w|V / |Vw|^, 
-'Q2 -'Q2 -'Q 

to show inequality (54), it is sufficient to prove the following inequalities 



Vw| dy < c ( I |e(w)|^ dy ) , (56) 

Q 



/ |Vwi| dy < c I / |V-wi| dy J . 

JQ2 \JQ2 J 



(57) 



Inequality (56) directly follows from Lemma 4.5. (iii) and Lemma 4.6. 

Let us now show that inequality (57) holds: let Wn G C^{Q2) be such that 
w„ —7- wi strongly in H^. Then, by integration by parts and Lemma 4.4 

Vw„Vwi dy = -(w„, Awi) = -(w„,Vv9) 
99V • Wn dy 



(58) 

^ IIv'IIl2(Q2) IIv- w„||^2(q^) 



< c ||v • w„|Il2(q^) ( W'^vW H-HQ2) + 



if dy 
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\Aw\ 



jj-x = sup 

MhHQ2)=^ 



Vw • Vu 



Q2 
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< sup ||Vu;||^.(Q^) l|Vn||i.(Q^) 

ueH^{Q2) 
ueH},{Q2) 

r"llifi(Q2)~"'" 
which imphes A defines a bounded Unear operator from H^ to H^^. Therefore 

Without loss of generahty, because adding a constant to (p does not affect Vc/?, we 
can choose Jq c^ dy = 0. Hence (58) becomes 



Vtt;„Vtt;i dy < c ( / |V • tf„| dy 

'Q-2 



\Vwi\^ dy 



Passing to the hmit n — )■ oo gives the desired result: 



Q2 



|Vtfi| dy<c\ I |V-ti;i| dy 



\Vwi\ dy , 



Q2 



implying (57). 



n 



4.3. On the technical lemmata utilised in the proof of Lemma 4-3 

Lemma 4.4: (see also Ladyzhenskaya [23].) There exists a constant c > such 
that for all u G L?'{Q2) 



I^IIl^(Q2) - ^ { II^"IIh-i(Q2) + 



u dy 



Proof: Assume the contrary. Then there exists {n„} in L^ such that ||iin|lL2(n \ 
1 for all n € N, and 



1 - \\un\\L2(Q^) > n [ ||Vn„||^_i(Q^-) + 



Un dy 



(59) 



Since L'^{Q2) ^^ H~^{Q2) is a compact embedding, {un} has a convergent subse- 
quence in H~^ to a limit uq say. After passing to the necessary subsequence we 
have 



VjTi 



uo in H 
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By Lions Lemma, c.f. [24, Theorem 3.2, Remark 3.1, p. Ill], 

\\Un - Um\\L^ < c{\\\7Un\\fj-i + ||VUm,||^_i + ||n„ - Um.\\H-i) n,m G N, 

which via (59) imphes {tt„} is a Cauchy sequence in L^ and, therefore, by the 
completeness of L^, converges to a Umit u say. Since L^ is embedded into H~^ we 
have 

ll'Un — "ullji^-i < c \\un — u\\^2 — > as n — > oo, 

which imphes u = uq, and also Jq Un dy — > Jq uq dy. Furthermore, Vn„ — > Vuq 
in H^^ since: for i = 1, • • • , n, 

I {Un,i -Uo^i,v) I = \-{Un -Uo,V^i)\ 

< C ||f,j||^2 \\Un - Uo\\^2 Vv G Hq{Q2). 



From (59) we see ||Vti„||j:^_i — > 0, Jq Un dy 



as n — > oo. Therefore 



||Vno||^-i =0, jQ^uody = 0, 

which implies uq = 0. Hence a contradiction as, from the construction of {un}, 
||no||i2 = 1. □ 

Lemma 4.5: For u G H],{Q) there exists u £ H],{Q) and a constant c > 
independent of u such that 

(i) u = u in Qi, 
(a) An = in Q2, 

(Hi) ||Vu||^2(Q) <c||u||^^(Q^). 

We shall call u the harmonic extension of u. 

Proof: For fixed u G H\{Q), Sobolev Extension theorem, c.f. e.g. [22], says that 
there exists an extension operator E : H\{Qi) — ;• H\{Q) such that 

||^u||//i(Q) <c||u||^i(Q^), (60) 

for some constant c > independent of u. Denote by tt G H^{Q) the solution to 

—Au = in (^2, u = u on dQ2, 

extended by u into Qi; u satisfies (i) &: (ii). Since u minimises the functional 

F{u) = / \Vu\^dy 

on {u + Hq{Q2)} we have, in particular, 

/ |Vn|2 dy< f |V(-Eu)p dy. 
JQ2 JQ2 

Inequality (iii) follows from (60). D 
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Lemma 4.6: Let u G [^i(<5i)] ; ^-e. u G [H^{Qi)] and Q-periodic. Then there 
exists a constant c > independent of u, such that 



H\m(Qi) ^ c 



u dy 
Qi 




Proof: As in the proof of Lemma 4.4, using the standard Korn's inequahty in the 
place of Lions lemma c.f. [24], we construct a subsequence ti„ that converges weakly 
in H^ and strongly in L^ to some uq, ||no||^2 = 1 with f^ uq = and e(no) = 0. 
Now, e(uo) = implies uq is a rigid body motion but, since uq is periodic, uq is not 
a rotation. Therefore, uq is constant, in particular uq = and the contradiction 
follows. n 



5. On the spectral convergence 

In this section we shall prove the spectral compactness result. Theorem 2.3. First, 
we note the validity of the result is subject to the following geometric alteration: 
the inclusions intersecting or touching the boundary are "excluded", i.e. we re- 
define C"'(x) to equal C^ (-) on such inclusions. Next, we give a precise definition 
as to the meaning of A^ and A^. 

Denote by H the closure of [H^in)f x [L'^{n; H^{Q2))f in L'^{n x Q) with inner 
product 

[u,v)^= I I u{x,y) ■v{x,y) dydx. 
Jn JQ 

Theorem 2.1 implies that 

A^{u,v) ■■= [ [ C^°'"VxuVxV dydx + [ [ a^^'^VyU ■ VyV dydx, 

defines a bilinear form on the Hilbert space H. The bilinear form A^ is, clearly, 
non-negative, and is closed on H. Therefore, it is well known, see e.g. [25] , that A^ 
defines a non- negative self-adjoint operator A^ , called the Friedrichs extension, by 

A^{u, v) = (^%, v)h, Vu, V G P(A°), 

where the domain 'D{A^) is a dense subset of H. We call A^ the homogenised limit 
operator. Similarly, we denote by A'^ the Friedrichs extension associated to the 
bilinear form (22). 



5.1. Convergence of spectra 

Lemma 5.1: The spectrum of A^ , a{A^), converges in the sense of Hausdorff to 
the spectrum of A^ , a{A^). That is 

(i) For every A G <y{A^) there exists X^ G cr(A^) such that A^ — )■ A. 
(a) If there exists A^ G o'{A^) such that A"^ — )■ A, then A G a{A^). 

According to Lemma 3.6, property (i) is implied by the strong two-scale resolvent 
convergence of A^ to A, which was proven in Theorem 2.1. We shall prove property 
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(ii) by arguments that are conceptually similar to [10]. Assuming A^ — t- A'^, let u'^ 
be the corresponding normalised eigenfunctions of A*", i.e. 

A'u' = \'u', lk'llL2(n) = 1- (61) 

Since the sequence u'^ is bounded, u'^ -^ u^ for some u^ £ L^(il x Q). By Theorem 
(2.1), we can pass to the limit e — )■ and find 

A%^ = A%°, 

To assure A*^ is in the spectrum of A^ is to show u^ is not identically zero. 
The quadratic form 

Biu, v) := / Vyuiy) ■ Vyv{y) dy (62) 

JQ2 

on the domain H := {v (i Hq{Q2) : V • w = 0} defines a self-adjoint operator B. 
It is well known, cf. [20], the operator B has a compact resolvent and therefore 
a discrete spectrum of eigenvalues going to infinity. Furthermore, by Lemma 5.5 
a{B) = {fin : n e N} C o-{A^). To prove u° / it is sufficient to show the following 
strong two-scale compactness result 

Lemma 5.2: Suppose that 

A. U — A U , \\U ||i2(f^) — i. 

Let A^ — 7- A ^ ^{B)- Then Ue has a strongly two-scale convergent subsequence. 

Indeed, if Lemma 5.2 holds then, by Definition (3.4), ||^i^||2,2(Q) — > II^^IL^mxO) 
which implies u^ ^ 0. 

Proof: Let J := {j G Z'^ : e{Q + j) fl il}, and u'^\j to be the restriction of u^ 
to j. By the geometric alteration made at the beginning of Section 5 we can use 
Lemma 5.4 below, to construct a v!^\j for all j G J ^ and therefore construct a 
v!" G Hl{n) such that u= = uMn $7f ; V • u^ = V • n'' in 0.% and Au'' = V(/)^ 
for some 0^ G Lp'{Q!2)- A straightforward rescaling of Lemma 5.4 (iii), inequalities 
(23)- (25) and jjw'^ 11^2(^2) = 1 show ||'u'^||//i(q) < C. Therefore, a subsequence of if 
convergences strongly in L^(r2) to some u. To prove the result it remains to show 
the difference v^ = u^ — u^ strongly two-scale converges. 

By construction w^ G [Hq{^2)Y7 ^y • v'^ = and since u"^ solves 

/ C(^)Vw^ • V0 dx + e^ f C7(o)Vw" ■V(pdx= [ X'u'(p dx V0 G C^i^), 
Jn Jq.1 Jo. 

we see that f ^ solves 

£2 / sjv^ . V(/. dx = / X^iv" + u^)<j) dx, V(/> G C^{9.1) such that V • </) = 0. 

(63) 

Let us consider the following variational Stokes problem: Find w^ G [i7Q(J72)]'^, 
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Vy • ti;"^ = such that 

e^ / Vw' ■V(f>dx= I X^w'cp + f(t) dx, V(/> G [C^{^l)Y such that V • (/> = 0, 

(64) 

for a given /^. It remains to show the fohowing result: If /^ ^ / then w^ -^ w 
where w G [L"^ {f!,; Hq)Y solves 



Vyw{x,y) ■ip{x)Vy(j){y) dydx = / Xw{x,y)ip{x)(l){y) dydx 





+ / /(x,y)V^(x),^(y)dydx, y,p e C^{n),4> e [C^{Q2)]'',Vy ■ ^ = 0. (65) 
Indeed, if this is true then 



hm / /^u^ dx = lim / e^Vw'^Vv'^ — X^w'^v'^ dx = hm / X^vfw^ dx 

Xuw dydx = / / Vyf • VyW — Xvw dydx 
JnJQ2 



fv dydx, 



where the first two equalities come from choosing test functions (;^ = u"^ in (64) and 
(p = w'^ in (63) respectively; the third equality uses the fact that u^ ^ u and finally 
we use symmetry of limiting problem (65). By choosing /^ = v^ we have shown, in 
particular, that 

lim / {v^{x)) dx = / V {x,y) dydx. 
^-^0 Jo. Jn Jq2 

Hence v^ ^ v. 

To show (65) we note that since the domain Ql consists of disjoint balls, the 
spectrum of B^, the variational Stokes operator defined on the physical domain 
r^l by (64), coincides with the spectrum of B on a single isolated inclusion as 
defined via (62) (by change of variables in (64)). Therefore, since {A^} is a bounded 
sequence and A"^ ^ '^(B) for small enough e, 

11^ ^ 11^^ - p(A^a(.4o)) - 



where p is the distance function. Hence, since w^ solves (64), 
W^^Wlhu) = \\iB'->^T'rh^n) <C\\f\\LHn) 



Furthermore 



/ \eVw'\'^ dx= [ X'iw'f + fw' dx < C. 
Jq.% Jo. 
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By standard two-scale convergence arguments, we see that w^ -^ w{x,y) G 
[L'^ {Q; H)]'^ , eVw'^ -^ Vyw{x,y). Passing to the two-scale limit in (64) gives (65). 

D 

Remark 1 : By arguments similar to those in [11] one can appropriately modify 
the above proof to show Lemma 5.1 for the case Q, = R'^. 

Lemma 5.1 tells us in particular that for small enough e the bottom of the 
spectrum of A'^ is a positive distance away from zero. Taking into consideration 
the Spectral Theory for self-adjoint operators we can see that a consequence of this 
is the following Poincare-type inequality: 

Lemma 5.3: There exists a constant C > 0, independent of e, such that for all 
u G [Hq{^)] and small enough e 



|^|lz,2(f^) < C 



2 



L^{ni) 



cW(x/e)Vu +e^ C(°)(x/e)Vn . (66) 



2 



L2{Q|) 



Proof: Since the bilinear form 

A^{u,v):= I (c^^\x/e)Vu{x)-Vv{x) + e^C^^\x/e)Vu{x)-Vv{x)\ dx 

•J i,l 

defines a non-negative quadratic form on the Hilbert space [-ffg(O)] , there exists 
a corresponding self-adjoint operator A^ such that on its domain T){A^) which is a 
a dense subset ViA") of [iJgY 

A^'iu.v) = (yl^ii,t;)i2(^). 

By the classical Rayleigh variational principle, 

A'^{u,u) = {A^u^u) > Xq{u,u), 

where Aq is the smallest eigenvalue of A^. By Lemma 5.1, Aq — )■ A'^ as e — )■ where 
Aq > is the smallest eigenvalue of the limit spectrum a{A^). Therefore, 

A'^{u,u) > c{u,u), 

for some c > and for all u G T){A^). This above inequality holds on [Hq] by the 
fact V{A^) is dense in [H^f. Hence, (66) holds. D 

Lemma 5.4: There exists a constant c> 0, such that for all u G [H^{Q)Y there 
exists u G [H^{Q)Y', such that 

(i) u = u in Qi, 
(ii) V • ii = V • n in Q2, 

(Hi) ||n||^i(Q) < c (||n||^.(Q^) + ||e(n)||^.(Q^) + ||V • uWl^^Q^^ , 

(iv) Au = \7(/) in Q2 for some (p G Lp'{Q2)- 

Proof: Introducing the space 

U = {ve [H^{Q)Y ■■ v{y) = k + X2{v)w{v) for some A: G R^ u; G [i^o (^2)]^ 

VyW = Q}, (67) 
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we introduce U , the orthogonal complement to U with respect to the following 
equivalent H^ {Q)-noiiai 

\Mh{Q) = IMl^Qi) + I|V^IIl2(q) • 

For fixed u G [^i(Q)]^; u = ui + U2 for ui ^ U and U2 G t/"*". Define f := n2 + A:, 
where k is the value of ui in Qi- Clearly t; satisfies (i) and (ii), as for y £ Qi, 

v{y) = U2{y) + k = U2{y) + ui{y) = u{y), 
and, in Q2, 

Vy ■ V = Vy ■ U2 =Vy ■ U2 + Vy ■ Ul =Vy ■ U. 

It remains to show (iii). This can be seen by following the proof of Lemma 4.3 
and using, where appropriate, the standard Korn's inequality instead of Lemma 
4.6. D 

5.2. Spectrum of the two-scale homogenised limit operator 

Let us study the eigenvalues of the homogenised limit operator A^ determined by 
Theorem 2.1. As outlined in Section 2 this requires studying the system 

-V • (&'°"'Vu{x)] = Xu{x) + X<v> (x) in f], (68) 

-Ayv{x,y) = Xv{x,y) +Vyp{x,y) in Q2 
Vyv{x,y) = inQs (69) 

v{x,y) = on dQ2, 

for some unknown p £ L'^{fl; H^{Q2)). 

For the elasticity problem with elasticity tensor (7'^°™, it is well known that the 
spectrum of the Dirichlet problem in Q is discrete and consists of eigenvalues A^, 
n = 1, 2, . . . , 00, such that 

< Af < Af < Af < . . . Af ^ 00 

with associated eigenfunctions u„ G [Hq{^)]'^ such that 

-V • (c^°'"Vu„(x)) = X^unix) in n. 

By setting {v,p) = (0,0) we see that A„ are in the spectrum of A^ with corre- 
sponding eigenfunction u^{x,y) = Un{x). 

As is also well known, for the Stokes spectral problem (69), the spectrum is also 
discrete and consists of Hm, m = 1, 2, . . ., such that 

< ^1 < /U2 < /is < • • • /Um, -^ 00 
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and Vm G [i?o((52)]'^, Pm G H^{Q2) such that 

-AyVmiy) = tJ-mVmiy) + "^yPmiv) in Q2 

'^y-Vm{y)=0 in (52- 

If, for some m, (vm) = then Hm is clearly in the spectrum of A^ with correspond- 
ing eigenfunction u^{x,y) = Vm{y)- For the eigenvalues ^m whose corresponding 
eigenfunctions have non-zero mean, i.e. {vm) / 0, assuming /x^ / ^n fo'^ ^li " l^t 
Um{x) G [i?d(ri)]'^ be the solution of 

-V • f C'^°"'VUm,(x) j = ^irnUm{x) + fJ-m < Vm > {x) in O. 

Then /i^ is in the spectrum of A^ with corresponding eigenfunction u^{x,y) = 
Um{x) + Vmix, y). If fim = Ki for some n we have already shown above that Hm he 
in the spectrum of A^. 

It remains to show that this exhausts all possible eigenvalues of the spectrum of 
A^. That is the following result holds. 

Lemma 5.5 Spectrum of the limit operator The spectrum of the homogenised 
lim,it operator A^, a{A^), has the following representation: 

a{A^) = {Xn I n G N} U {nm \ m G N}. 

Here Xn satisfies, for some non-trivial Un G [Hq{^)Y' , 

-V • {c^°'^VUn{x)\ = XnUn{x) in fl, 

and fim satisfies, for some non-trivial Vm G [-^^0(^2)]'^; Pm G H^{Q2), 

-AyVmiy) = l^mVm{y) + "^ yPm{y) in Q2 

^y • Vrn{y) =0 in Q2- 

Proof: For A„, ^m given in Lemma 5.5 we have shown that 

{A„ I n G N} U {nm \m G N} C cf{A^). 

To show the reverse inclusion it is sufficient to show that if A 7^ A„, A 7^ ^m; 
Vn, Vm then for a given f{x,y) G [L'^{Q, x Q)Y there exists a unique solution 
u = uq{x) + v{x,y), uq G [Hq{Q)]'^, v G [L'^{Q,;Hq{Q2))Y, continuously depending 
on / to 

-V^ • (^C^°'^VMx)) - Mx) = X<v> {x)+ <f> (x), (70) 



and 



-Ayv{x, y) - Xv{x, y) = f{x, y) + Vyp{x, y), 
-Vy ■v{x,y) = 0, 
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for souiep{x, y) G [L?'{Q.; H^{Q2))] ■ It is well known that if A 7^ ix^m Vm, then there 
exists a unique v to problem (71). Furthermore, g := \ < v > + < f >G [L^(r2)] 
and it is well known if A ^ An, Vn, there exists a unique uq G [/7Q(r2)] such that 



-V, • (c^°"^V,n(x)) - \u{x) = g{x). 



The above construction ensures, by the boundedness of the appropriate inverse 
operators, the continuity of the solution in /. D 
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